Properties of low— lying states in a diffusive quantum dot and Fock-space localization 
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Motivated by an experiment by Sivan et al. (Europhys. 
Lett. 25, 605 (1994)) and by subsequent theoretical work on 
localization in Fock space, we study numerically a hierarchi- 
cal model for a finite many-body system of Fermions moving 
in a disordered potential and coupled by a two-body interac- 
tion. We focus attention on the low-lying states close to the 
Fermi energy. Both the spreading width and the participation 
number depend smoothly on excitation energy. This behav- 
ior is in keeping with naive expectations and does not display 
Anderson localization. We show that the model reproduces 
essential features of the experiment by Sivan et al. 

PACS numbers: 72.15.Rn, 73.23.-b 



The measurement of the quasiparticle spectrum of a 
diffusive quantum dot via its tunneling conductance by 
Sivan et al. Q] in 1994 has caused considerable theoreti- 
cal activity. The experimental spectrum displayed a few 
narrow peaks in the vicinity of the Fermi energy, followed 
by a quasi-continuum Q . The number of discrete peaks 
was found to be of the order of the dimcnsionless con- 
ductance g of the dot. The results have given rise to a 
debate on the validity of Fermi liquid theory for diffusive 
quantum dots. In this theory, the low-lying excitations of 
a system of interacting Fermions are described as quasi- 
particles, i.e., free Fermions with a renormalized mass 
and a finite life time r, visible as peaks of width Ti/t in 
the spectrum. Are these predictions consistent with the 
very limited number of peaks visible in the spectrum of 
Ref. [El, or does disorder invalidate Fermi liquid theory? 

Altshuler et al. || asked (and answered) this question 
in a precise way. These authors related the many-body 
problem with electron-electron interaction in a diffusive 
quantum dot to that of single-electron Anderson local- 
ization in real space. Many-body Fock states are intro- 
duced as Slater determinants of eigenstates of the single- 
particle Hamiltonian containing kinetic energy and dis- 
order potential. The distance in Fock space between 
two such states is defined as twice the minimal number 
of electrons which have to be moved from one single 
particle state to another in order to get from one state 
to the other. Matrix elements of a two-body interac- 
tion between two states differ from zero only if the dis- 
tance between the two states is < 4. After introducing 
the Fermi energy and a particle-hole representation, the 
states are grouped into classes, each class being defined 
by the number of particle-hole (p — h) excitations. The 
distance between classes is defined as the minimum dis- 



tance between any state in one class and any state in 
the other. In the work of Ref. ||, only couplings with 
distance 2 were considered. Moreover, only terms that 
increase the complexity of the states were kept, i.e., cou- 
plings to all states in the same and to all but one state in 
the next-lower class were neglected. With these assump- 
tions, the Fock-space problem could be mapped onto a 
tight-binding model on the infinite Bethe lattice. 

For g 1, the existence of three regimes separated 
by two characteristic energies was shown. The Ander- 
son transition [Q between localized and extended states 
occurs at the energy E** ~ Ay/gjhig, while E* ~ A^fg 
defines an effective energy above which the many-particle 
states are completely mixed. Here, A is the average 
single-particle level spacing. These results were corrobo- 
rated by a calculation using supersymmetry 0] and later 
also discussed in relation to the problem of few interact- 
ing particles in a random potential to the two-body 
random interaction model 0, and to the level statistics 
of excited many-body states || . In Ref. , doubts were 
voiced on the claim for a delocalization transition. In 
Ref. H, finite-size effects not considered in the calcu- 
lation using the Bethe lattice were estimated to yield 
E c \y ~ Ag 2/>3 for the value of the energy at which the 
states become completely mixed. 

The work of Refs. uses a number of approxima- 
tions which are needed to establish the connection with 
the Bethe lattice. In addition, it is assumed that the 
Bethe lattice has infinite length. These approximations 
are presumably valid at sufficiently high excitation en- 
ergy E. On the other hand, the data of Ref. refer to 
the immediate vicinity of the Fermi energy. For the inter- 
pretation of these data, it is important to know whether 
the results of the high-energy approximations apply. For 
instance, at low excitation energy and with e = E/A, the 
number of accessible classes has a strong cutoff ~ yfe. 

To answer this question, we investigate in this Letter 
the low-lying states of the many-body problem with dis- 
order and interaction in the framework of a model which 
is more realistic than but retains the spirit of Refs. [B,§J- 
We avoid the approximations made in these papers which 
are needed to obtain the structure of the Bethe lattice. 
We pay the price that we cannot use analytical approxi- 
mations and must rely on numerical simulation. 

The Hamiltonian for spinless electrons (Fermions) has 
the form H = Hq + V, with Ho = J2k £ fc a l a fc the un- 
perturbed Hamiltonian, a sum of single-particle opera- 
tors, and V the two-body interaction. The symbol 



1 



denotes the single-particle eigenvalues, and at generates 
the single-particle states \k). The single-particle Hamil- 
tonian contains a random potential. Therefore, the states 
\k) and energies have stochastic properties. In an 
energy interval of length <?A, these properties coincide 
with those of the Gaussian orthogonal ensemble (GOE) 
of random matrices. We use the classification scheme de- 
fined above and consider classes m of Fock states with 

m particles and (m— 1) holes. Here, m = 1,2, We 

denote the states by |m, i), where i is an index run- 
ning over the states of the class. The associated ener- 
gies, given by sums of the e^'s, are denoted by E myi , so 
that Ho\m,i) = E mt i\m,i) . To implement this model, we 
take the energies eu after unfolding from the center of 
the Wigner semicircle for the GOE. We choose the Fermi 
energy equal to zero. 

The unperturbed mean level density p^(e) in class m 
is given by || p° m {e) — A _1 e 2m /[(w + l)!™!(2rn)!], where 
e = E/A and E is the excitation energy. With s fixed, p® m 
grows strongly with m until it suddenly drops to almost 
zero at m ~ -s/e. Hence, at any value of e only a limited 
number of classes contributes to the total level density 
of the unperturbed system (V — 0). In the vicinity of 
the Fermi energy, this number is one or two. This fact is 
important for the experiment of Ref. Indeed, if the 
localization length is larger than two, then localization in 
Fock space can have no bearing on the spectrum in the 
vicinity of the Fermi surface. 

The interaction operator V mixes the states |m, i). We 
suppose that the diagonal part of the interaction is in- 
cluded in Hq by use of the Hartree-Fock method, without 
affecting the statistical properties of either the energies 
E m i or of the states \m,i). The matrix elements of V 
between two different Fock states \m, i) and \n,j) vanish 
unless |m — n\ = or 1 and unless the Fock distance 
of both states is < 4. (The case \m — n\ = 2 requires 
creation or destruction of two particle-hole pairs out of 
the Fermi sea. Such processes contribute to the unlinked 
diagrams of perturbation theory and are not considered 
here. Calculations including couplings with \m— n\ = 2 
yielded results which did not differ significantly.) The 
non-vanishing matrix elements of V are assumed PJlO| 
to have a Gaussian distribution centered at zero with 
variance V 2 = (A/g) 2 . The mixing of the states \m, i) 
depends on their spacing and on the strength of V, i.e., 
on the value of g. In our model, the strength of V does 
not depend on excitation energy. This fact will cause 
stronger mixing of the higher-lying states with smaller 
spacings. The effect is compensated because out of the 
larger number of states, a decreasing fraction couples to a 
given one. We allow for couplings between states in the 
same class. This introduces terms which are explicitly 
excluded in the case of the Bethe lattice. 

In our numerical work, the energy scale is defined by 
A. Our only parameter is the dimensionless conductance 
g. The dimension of the matrices was limited by a cutoff 
E cuto e^> A: All Fock states with energies E m ,i > E cu to& 
were omitted. For each realization of H, E cuto g was cho- 



sen in such a way that the matrix dimension was 1000. 
This corresponds roughly to -E cu toff ~ 18 A. Since our 
single-particle energies are drawn from the center of a 
GOE distribution, consistency requires that the Thou- 
less energy gA is at least as large as the energy interval 
considered. The eigenstates |a) and eigenvalues E a of the 
full problem were obtained by diagonalization. We have 
checked that the results presented below are independent 
of E cut off ■ We now present two statistical measures suit- 
able for a test of localization in Fock space. 
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FIG. 1. Quasiparticle line shapes for the lp— Oh states 
|1, 2), 1, 6), and 1, 11) (full curves). The dotted curve is the 
quasiparticle spectrum of all lp — Oh states. The inset shows 
the spreading width F versus E for g = 5 (circles), g — 10 
(squares), and g = 15 (crosses) together with the Golden Rule 
predictions Fgr(-B) (full curves). 



Quasiparticle spectrum. Adding an electron to a quan- 
tum dot with a filled Fermi sea creates a lp — Oh Fock 
state \l,i). The interaction V spreads this state over a 
set of eigenstates \a). The spectral shape of the resulting 
quasiparticle peak, also referred to as the local density of 
states, is given by 



Pi,i(£) = (£l<«lM>| a S(E-E a ) 



(1) 



The brackets denote the ensemble average. In order to 
obtain meaningful plots, we have used only a single real- 
ization of the single-particle energies but 100 realizations 
of the interaction matrix elements. To reduce the den- 
sity of the resulting data set, we have averaged our results 
over small energy windows. This yields the dots shown 
in Fig. |l| for three selected lp—Oh states at g = 10. The 
fits with the Lorentzian 
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2n (E-e hi ) 2 +T 2 /4 



(2) 



(full curves) fix the centers of the quasiparticle peaks 
and the spreading widths r^j. Adding the Lorentzian 
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line shapes of all Xp — Oh states yields the quasiparticle 
spectrum shown as dotted line. 

Spreading width. Plotting Tij for all states ver- 
sus the state energies e±^, we obtain the spreading width 
T(E) as a function of excitation energy. For several values 
of g, this dependence is displayed in the inset of Fig. [j]. 
We compare our results with the prediction of the Golden 
Rule. We study the mixing of lp—Oh states. Each one of 
these is directly coupled to both all other lp— Oh states 
and all 2p—\h states. The lp—Oh states have mean spac- 
ing A. At the coupling strengths (g>5) shown in Fig. []], 
these states are barely mixed with each other. There- 
fore we keep only the density of 2p—lh states, P2(E), 
as the appropriate quantity to use in the Golden Rule 
expression 

T GR (E) = 2tt (A/g) 2 p 2 (E). (3) 

We point out that this relation is valid beyond perturba- 
tion theory if the density of final states is taken to be the 
exact (rather than the unperturbed) level density. The 
exact density in class m is defined as 

Pm(E) = /^|(a|m,i)| 2 6{E-E a )\ . (4) 

\ i,a I 

The brackets again denote the ensemble average which 
in practice we perform over 100 realizations of the full 
Hamiltonian H. The three full lines in the inset of Fig. |l| 
show TgrX-E) for the three values of g. With decreas- 
ing g, i.e., increasing interaction V, Tq R (E) increases 
while the slope of Tgb.(E) decreases. These facts re- 
flect both the 1/g 2 dependence of V 2 and the behav- 
ior of P2(E). For g = 10, we show in Fig. ^ Pi(E), 
P2(E), Pt{E) — J2 m Pm(E), and the unperturbed den- 
sity of states p% (E) . The latter differs substantially from 
P2{E): The interaction causes the density to become 
wider, thereby reducing its slope. This accounts for the 
behavior of r GR (^) in the inset of Fig. |. In Refs. §,§, 
it was emphasized that localization in Fock space invali- 
dates the Golden Rule. The quantitative agreement be- 
tween our numerical results and the expression (^) down 
to very small energies shown in Fig. [l] is, thus, a very 
strong argument against the occurrence of localization in 
the energy range investigated in this paper. The weak 
mixing of the low-lying states is entirely explained by 
the low density of states right above the Fermi energy 
and is not due to additional restrictions in Fock space. 

Qualitative features of the quasiparticle spectrum in 
Fig. [l] are in good agreement with the results of the 
experiment by Sivan et at 0J^]: A few nearly discrete 
quasiparticle peaks with T < A occur right above the 
Fermi energy. The width T grows with excitation energy 
until r « A at E w gA where the spectrum becomes 
quasi-continuous. For an accurate simulation of Sivan's 
experiment, however, the lp — Oh excitations should be 
built upon the ground state of interacting particles rather 
than upon the filled Fermi sea. Calculations using such 



an interacting ground state reproduced the results shown 
in Figs. [I] and except for g — 5. Here the spreading 
widths were at the upper end of the statistical fluctua- 
tions shown in Fig. H, while the Golden Rule underesti- 
mates r by about a factor two. The clear hierarchical 
structure needed to study localization is lost when we 
start from an interacting ground state. Therefore, we 
kept the simple model outlined above. 
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FIG. 2. Densities pi{E), p 2 (E), and p t (E) for g = 10 as 
defined in the text. For comparison we also show p%(E). 

Participation number. The average participation num- 
ber often serves as a measure of localization. The parti- 
cipation number for the Fock state \m,i) is defined by 

a 

If the state \m,i) is localized, R m ,i ~ 1, while R m ,i in- 
creases monotonically with increasing mixing. In the 
thermodynamic limit, R m ^ is unbounded. In the case 
of finite matrix dimension N and the GOE, there exists 
an upper bound i? ma x = (1/3) TV for i?„ v ;. In Fig. || we 
show the average participation number R\(E) of class 
1. This quantity is obtained by plotting the values of 
Ri i for all available i versus the corresponding energies 
Ei i f° r a number of realizations and averaging the re- 
sult over energy. For weak interaction (here g = 15) and 
low excitation energies, the Fock states \ l,i) are localized 
because the spreading widths of the quasiparticle states 
are small compared to the mean level spacing. With in- 
creasing V (g = 5) the mixing becomes larger even at 
low excitation energies. We only find smooth transitions 
to mixed states, both with increasing excitation energy 
E and with increasing interaction V. From the absence 
of any discontinuity in R\{E), we conclude that there is 
no evidence for a localization transition in Fock space. 
For g = 15, the derealization thresholds predicted in 
Refs. §,§ have the values E** « 2.34A and E* « 3.87A. 
We find R\(E) < 2 for E < 5 which shows that at these 
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energies, the mixing of Fock states is rather weak. More 
generally, an estimate of the number of eigenstates \a) 
contributing to a Fock state at energy E and, thus, 
of the average participation number R\(E) is given by 
Rf R (E) = T GR (E) p t (E) (see Eq. (§) and Fig. |). We 
find this relation to be correct for sufficiently high ener- 
gies E (see Fig. ||) while Rf R (E) overestimates Ri(E) at 
low energies since by definition, R\ (E) > 1 whereas the 
level density becomes very small. For a localization tran- 
sition, we would expect Ri(E) to lie below the Golden 
Rule estimate. This is not the case. 
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FIG. 3. Participation number Ri (E) of the lp-Oft states for 
g — 15 (full circles) and g = 5 (empty circles). The estimate 
Rf R (E) for g = 15 is plotted as a full curve. In addition the 
results for applying a loopless Bethe lattice model (crosses) 
and considering couplings over Fock distance 2 only (stars) 
are shown (both g = 15). 

Sensitivity to coupling scheme. How do our results 
change with coupling scheme, i.e., with the omission of 
those interaction terms which spoil the analogy with the 
Bethe lattice? In the coupling scheme applied so far we 
allowed for couplings between states of Fock distances 2 
or 4, in keeping with the assumption of a two-body inter- 
action. The stars in Fig. [|show the average participation 
number for Ip — Oh Fock states when only couplings be- 
tween states of Fock distance 2 are taken into account. 
This restriction means that only one particle is allowed 
to change its single-particle state in an interaction and 
no other particle can compensate the energy difference 
implied by this transition. In effect couplings between 
states close in energy are neglected. 

In order to test the Bethe lattice assumption, i.e., a 
Fock space topology without any loops, we first assumed 
that every Fock state \m, i) with to> 2 is coupled to only 
one state in class to— 1 . This state was taken to be the one 
closest in energy to |m, i). Couplings within each class 
were neglected. Couplings with Fock-space distance 4 
were taken into account. This model yielded a set of dis- 
connected trees each starting from one lp — Oh state. As 



each Fock state belongs only to exactly one tree, the trees 
were all of rather small size. The resulting participation 
numbers were much smaller than those of the full calcu- 
lation, especially at high energies. This is unsatisfactory. 
We then considered a second model. We coupled each 
lp — Oft, state to all 2p—lh states and kept the couplings 
between higher classes as previously. Now all Fock states 
except the lp— Oh states participate in each tree. The re- 
sulting average participation numbers (crosses in Fig. |^) 
agree rather well with the full calculation. 

In conclusion, we presented a numerical study of 
the eigenvalues and cigcnfunctions of a finite Fermionic 
many-body problem with random single-particle ener- 
gies and a random two-body interaction. Our model 
avoids the simplifications of the Bethe lattice. Matrix 
diagonalization with a cutoff yielded results restricted to 
the vicinity of the Fermi surface and insensitive to the 
choice of the cutoff. We found that it is important to 
keep couplings between states with a distance 4 in Fock 
space. Moreover, we showed that for sensible values of 
the conductance <?, our model is able to reproduce es- 
sential features of the experiment of Ref. |ij . Calculated 
values of the spreading width and the participation num- 
ber indicate some degree of localization, depending on 
interaction strength and excitation energy. However, in 
contrast to analytical predictions based on a high-energy 
approximation and on the Bethe lattice (which claim a 
twofold localization transition) as well as those of Ref. [|| , 
we have only found a smooth transition from almost lo- 
calized to delocalized states. This transition can be fully 
understood in terms of the density of available lp— Oft. and 
2p— lft states and does not display Anderson localization. 
Our results show that the behavior of a finite Fermi sys- 
tem at low excitation energy and at zero temperature 
differs profoundly from the thermodynamic limit. 
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